NASA TM X- 68193 


NASA TECHNICAL 

MEMORANDUM 


N73- 19299 
NASA TM X- 68193 



APPLICATION OF THE PRINCIPLE OF CORRESPONDING 
STATES TO TWO-PHASE CHOKED FLOW 


by Robert C. Hendricks and Robert J. Simoneau 
Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PAPER proposed for presentation at 
Seventy- fourth National Meeting of the American 
Institute of Chemical Engineers 
New Orleans, Louisiana, March 11-15, 1973 



E-7340 


APPLICATION OF THE PRINCIPLE OF CORRESPONDING 
STATES TO TWO-PHASE CHOKED FLOW 


Robert C. Hendricks and Robert J. Simoneau 
Introduction 

The state of research in two-phase choked flow, while sub- 
stantially developed in recent years, still has significant unanswered 
questions. Hsu (ref. 1) presents a good review on this. The proper 
choice of fluid mechanic and thermodynamic models is not totally 
resolved (refs. 2 and 3). Henry and Fauske (ref. 3) make a strong 
argument that nonequilibrium effects must be considered; however, 
in doing so they are forced into empiricisms. As a result, designers 
are reluctant to work with the theories beyond the range of known data, 
and frequently demand data for the fluid they are using and in the range 
they are working. In the case of oxygen this is particularly distressing 
since testing is expensive and hazardous. 

With this in mind the present paper examines the possibility of 
applying the principle of corresponding states to two -phase choked 
flow. If this is valid, designers, using data from other fluids to 
back them up, could apply the theory of their choice to fluids where 
there is little or no data, such as oxygen. 

In ref. 4 it was noted that, while the isenthalpic Joule -Thomson 
coefficient (jj) is not a corresponding states parameter, a reducing 
parameter (ju, c T c /6P c ) could be defined and several figures of 
li R = jut/ju c as a function of T R = T/T c with P R = P/P c as a parameter 
were presented. Several fluids including methane, oxygen, and nitro- 
gen appear to form an average parametric plot which indicates that 
I a must nearly obey the principle of corresponding states as discussed 
in ref. 4. With this in mind, it was assumed that there could be 
several thermodynamic flow processes which nearly obey the principle 
and consequently could be placed in some reduced form. It was decided 
to examine whether two -phase choked flow could be one of them. 
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As a cautionary note, the authors wish to emphasize that in 
this paper they are not advocating any particular theoretical analysis 
for two -phase choked flow. In fact, the authors believe that when 
the final verdict is in, it will involve a combination of theories applicable 
in various regions. The purpose here is simply to evaluate corresponding 
states normalization. 

ANALYSIS 

While the arguments leading to this Conclusion are beyond the 
scope of this paper, it is well known that one-dimensional isentropic 
flow through a nozzle is related to the enthalpy difference. Since the 
enthalpy difference nearly obeys the corresponding states principle, 
it follows that nozzle flow should also nearly obey the principle. Such 
is the case even when fuel rich combustion is taking place as in a rocket 
engine, a well known result in rocketry. 


Critical Flow Rate 

The maximum flow rate (choked flow) of a gas through a nozzle 
has been shown by several authors (eg. ref. 1, 2, 5) to be a thermo- 
dynamic parameter dependent on the state points (P Q , T Q ) and (P t , T t ) 
which satisfy equation 1 
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In reference 6, the basic flow equation for both nonequilibrium 
and equilibrium models is given as 
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and the maximum flow rate was found by maximizing equation 2 
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By simply normalizing equation 2, the reduced maximum flow 

rate becomes: ^ ^ 

* „ 

■R (4) 


(2 — * 




If, however, the enthalpy form of the flow equation, equation (1), is 
used, the normalizing constants are somewhat different. In arguing 
that the enthalpy nearly obeys the principle, the reduced enthalpy was 
found to be 0 

H "RX ' "f ( T R , V ^ O, J ^ (5) 


where f and g are rather complex functions and f obeys the 
principle while g does not. Examination of equation (1) in light of 
the reduced enthalpy of equation (5) suggests the proper normalizing 
parameter to be 

<5 = 



Clearly for the near spherical molecules, methane, nitrogen, 

argon and oxygen, it makes little difference which form of G is 

* 

used; however equation 6 should be the preferred form for G . The 
critical constants used herein are from reference 8 and given as Table I. 


Critical Pressure Ratio 

The isentropic expansion of an ideal gas through a nozzle leads 
to the critical pressure ratio being ^ 


T? 


A. 


Jr' 



(7) 


Equation (7) is not reduceable by corresponding states since 
Cp,Cvand are not corresponding states variables. Consequently 
this would lead one to believe that the critical pressure ratio cannot 
be reduced, since the asymptotic value cannot be reduced. On the 
other hand, if the low density ( p — 0) background values, Cp Q and 
Cv Q are subtracted out then Cp-Cp Q and Cv-Cv Q can be nearly 
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reduced by corresponding states. This is discussed by Francis and 
Luckhurst (ref. 7). Such results suggest that ^ in equation (7) 
could be reduced in similar manner. The authors* efforts to subtract 
out this background and normalize ^ were not successful. However, 
it was noted that in the near critical region there existed a correspon- 
dence between the critical pressure ratios for nitrogen and methane 
for the same reduced pressure and temperature. If the principle is 
to work, it should work in the near critical where background effects become 
small. 

It was decided then to examine the reduced critical pressure 
ratio for both cases: i) the dense fluid region where the normalizing 
parameter is unity 




(dense fluid) (8) 


and ii) the gas region where the normalizing parameter is the perfect 
gas asymptote, eq. 7, 


Jk.\ „ 



(gas) 


(9) 


The asymptotic values (P t /P Q )pQ usec * herein are given as Table II. 
This is an ad hoc action which will be discussed further in the results. 


RESULTS 


Comparison with Isentropic Equilibrium Theory 

Equation 2 was solved assuming an isentropic equilibrium expansion 

and maximized to give the critical flow rate (G_ 0 „) and critical flow 

max 

pressure ratio (P t /P Q ) for the corresponding states fluids nitrogen, 
oxygen, argon, and methane. The thermodynamic properties were taken 
from reference 8. The choice of the equilibrium model was made 
because the calculation is very straightforward, because no empiricisms 
are involved and because it seems to work well in the thermodynamic 
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critical region, where the authors are especially concerned. A 
better choice of models may be appropriate in some regions, 
especially when P Q is near saturated liquid, but the equilibrium 
model will illustrate the corresponding states question. The 
normalization follows equations 6 and 8 or 9. 

Critical flow . - By systematically varying the stagnation 
pressure, holding stagnation temperature constant, one can produce 
a map of reduced critical flow rate as a function of reduced pressure 
with reduced temperature as a parameter. Such a map of four 
corresponding states fluids is given as figure 1. A closer examina- 
tion of fig. 1 indicates the calculated values for nitrogen and oxygen 
are the same; however the values for argon and methane tend to 
deviate, higher and lower respectively. Furthermore the departure 


from a single isotherm increases from 2% for 



1. 1. The substantial increase between 


T R < 1.0 to 10% at 
T R = 1.0 and T R = 1. 1 


is at this point unexplained. The departure for T R 2s 1.2 appears 
to be about 13%. This degree to which the four fluids collapse to a 


single isothermal curve both strengthens and weakens the corres- 


ponding states assumption. It is strenthened in that it works well 


for Tr~ 1.0. It is disturbing, however, that there is a shift 
between T R~ 1,0 and T r = 1.1, a region where the background 
(or noncorresponding states effects) should be small relative to the 


total and should subtract out. 


The reader should be reminded that, even though the flow 

normalization has not perfectly collapsed all the fluids to a single 

* 

reduced curve, the normalization by G has brought fluids with 
widely varying flow rates for the same reduced conditions very 
close together. For example at Pr 0 = 3. 0 and T R =1.1 the values 
of G r for argon and methane are 1. 34 and 1. 22 respectively, a 9. 3% 
spread; whereas the values of G are 12, 600 and 6, 200 gm/cm^sec 
respectively, a factor of 2 difference. 

The abnormal small ’’bumps” in these curves (such as at T D = 1.0, 

Jrt 

LOCPr <1.2) represent changes in the flow patterns which are best 
explained in terms of pressure ratio; therefore, discussion is deferred 
to the next section, where the ’’bumps" are more pronounced. 
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Critical flow pressure ratio . - The critical flow pressure 
ratio normalized by equation 9 as a function of reduced stagnation 
pressure with reduced stagnation temperature as a parameter is 
illustrated for argon, oxygen, nitrogen and methane in figure 2. 

Only the oxygen curves are plotted for values of T R > 1. 1 to avoid 
confusion, since the curves are getting close together and they are 
folding over. As with the critical flow rate the general observation 
can be made that representing critical flow pressure ratio as a 
function of reduced stagnation pressure and temperature tends to 
bring the four fluids together. While this is very encouraging, there 
is still a 10-15% spread in some areas, just as with the flow rate 
curves. Again the T R = 1.1 isotherm shows substantial deviation. 

It is interesting to note, however, that if Pj./P 0 ) R is plotted with 
the dense fluid normalization, equation 8, the deviation picture 
reverses as shown in figure 3. In the dense fluid region and near 
the thermodynamic critical region, (T R = 1.0, 1.0 JSPj^Si 3.0) the 
Pf/P^R curves for the four fluids collapse very nearly to a single 
curve as a function of reduced pressure and temperature. On the 
other hand, in the more gaseous regions (T R =1.4 or Pr q < 1.0) the 
spread appears again. Obviously one could contrive a more empirical 
normalization which would cover both regions, however, since the 
whole approach on pressure ratio is somewhat ad hoc the authors 
prefer not to do this and leave it to the discretion of the reader. 

As mentioned before, the "bumps” in these curves represent 
changes in flow patterns. For example, the peak in the T R = 1 isotherm 
corresponds to the first appearance of net quality in the nozzle. The 
inflection point near Pr q = 1.02 occurs when the stagnation entropy changes 
from liquid to vapor. The near horizontal region (P R < 1) represents 
combined two -phase and vapor flow (mist flow) and the final drop 
represents the transfer from two -phase to single phase flow throughout 
the nozzle. The horizontal line (at low pressure) represents single 
phase flow. 
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Comparison with Nitrogen and Methane Data 

Choked flow data were taken for fluid nitrogen and methane in 
a converging-diverging axisymmetric nozzle. The complete set of 
nitrogen data and description of the experiment are presented in 
reference 6. Subsequently methane data were taken in the same 
facility in order to make the corresponding states comparison pre- 
sented in this paper. From the data acquired for the two fluids it 
was possible to obtain four matching reduced isotherms for comparison. 
These were T^ = 0.81, 0.88, 0.96, and 1.035. The range of data 
obtainable was dictated by the test facility which was originally a 
liquid hydrogen heat transfer facility. The data are presented in 
Table III and shown, along with the equilibrium theory curves, on 
figures 4 and 5. 

A careful examination of figures 4 and 5 and T able in reveals that 
in general the nitrogen and methane experimental data reduce very well by 
corresponding states. If the data follow the equilibrium theory fairly 
well, such as along Tr = 1. 035, then both fluids do. If the data deviates 
from the equilibrium theory, such as along T^ = 0.96, then both fluids 
deviate — and in the same direction and by about the same amount. As 
stated earlier the authors do not wish to attempt to verify any parti- 
cular theory in this paper; however, some comments are in order. 

The data of this paper, and particularly that of reference 6, suggest 
that an isentropic equilibrium analysis will work well above T R ~ 1 '°* 

On the other hand, near saturated liquid stagnation conditions nonequili- 
brium effects are probably important and analyses, such as Henry and 
Fauske (ref. 3), will be required. The region in between is not clear 
yet and more experimental and analytical work is required. Substantial 
deviations in slope, such as the flagged data points along T R = 0. 88, 
cannot be correlated by either type theory at present, and may represent 
other factors such as two-dimensional effects in the nozzle. 
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CONCLUSIONS 

Two-phase critical flow and critical flow pressure ratio of 
simple fluids through a nozzle have been shown to nearly obey the 
principle of corresponding states. 

Normalizing the critical flow by 



nearly reduces the calculated equilibrium value of critical flow to a 
single isothermal curve for the fluids nitrogen, oxygen, methane 
and argon. 

For dense fluids and near the thermodynamic critical point 
the critical flow pressure ratio, P^/P 0 , is approximately the same 
for each fluid as a function of reduced stagnation pressure along lines 
of constant reduced stagnation temperature. In the less dense or 
gaseous region a normalization of critical flow pressure ratio by 

-f-) - f-H#r 

r ° PC L P ‘ 

seems to be required. 

In the region around T r = 1.1 the spread among the curves, 
which represents the degree of failure to follow corresponding states, 
is fairly large (greater than 10%). This is unexplained at present 
and requires further examination. 

The corresponding states normalization presented is also sup- 
ported in this paper by nitrogen and methane two-phase choked flow 
data. 

Since the reduced curves for nitrogen, oxygen and methane are 
so close, it is felt that oxygen two-phase choked flow can be predicted 
using corresponding states reduction and relying on nitrogen and methane 
data to establish the appropriate theory to use. Consequently only a 
limited number of experimental oxygen points are required to verify 
this conclusion. 
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NOMENCLATURE 


C 

Cp 

Cv 

C ?o * Cv o 

f,g 

G 

G* 

H 

P 

R 

T 

V 

X 

z 

r 

f 

Subscripts 

c 

g 

1 

max 

o 

PG 

R 

t 


Sonic velocity 

Specific heat at constant pressure 
Specific heat at constant volume 
Zero density specific heats 
Functions 
Mass flux 

Mass flux normalizing constant 

Enthalpy 

Pressure 

Gas constant 

Temperature 

Specific volume 

Quality 

Compressibility coefficient 
Specific heat ratio Cp/Cv 
Joule-Thomson coefficient 
Density 

Thermodynamic critical point 
Gas (vapor) 

Liquid 
Maximum 
Stagnation 
Perfect gas 

Reduced or normalized (see text) 
Throat 


cm/sec 
j/gm K 
j/gm K 
j/gm K 

gm/cm 2 sec 

gm/cm 2 sec 

j/gm 

MN/m 2 
j/gm K 
K 

cm 3 /gm 


K m 2 /MN 
gm/cm 3 



10 . 


FLUID 

P c 

MN/m 2 

T c 

K 

fc 

gm/cc 

Z 

c 

Z c 

G* 

gm/cm^sec 

Nitrogen 

3.417 

mm 

.3105 

.2937 

6010.4 

Oxygen 

5.083 

| 

.4325 

.2922 

8673.9 

Methane 

4. 627 

190.77 

.162 

.2889 

5093.7 

Argon 

4.865 

150.7 

.531 : 

.2921 

9404. 2 


Table I Critical Constants Used in the Reducing Parameters 


FLUID 

r 

p tl 

V 

PG 

Nitrogen 

1.4 

.5283 

Oxygen 

1.396 

.5290 

Methane 

1.304 

.5450 

Argon 

1.668 

. 4870 


Table II 


Asymptotic Values of the Critical Pressure Ratio 






















TABLE m NITROGEN AND METHANE CHOKED FLOW DATA 


11 


i-H E- O CO 

to O OS 03 


NOinnNOlOtDffiOlt- 
t- t- t- N 05 W CD- t— 00 H 


• 

tH 

CM* 


• 

• 

CO 

• 

E- 

* 

C- 

00 

• 

05 

oo 

oo 

05 

CO 

t- 


tH 

iH 

1— 1 

rH 









rH 

rH 

CO 

os 

eo 

rH 

oo 

o 

CO 


05 

rH 

e— 

l-H 

00 

CD 

00 

00 


CO 

os 

05 

CO 

CO 

CM 

E- 

CM 

e- 

E- 

CO 

00 

rH 

o 

in 

LO 

eo 

CO 


in 

in 

co’ 

•Cfi 

CM 

CM 

»H 

• 

o 

CM* 

• 

rH 














1 

rH 

6 


oo 

e- 

co 

in 

00 

oo 

CM 

00 

O 


05 

o 

05 

CO 




co 

t-H 

CD 

oo 

CM 

E- 

00 

CO 

CM 

rH 

CO 

in 

CM 

CM 



CM 

CM 

CM 

CO 

co 


in 

CD 

oo 

00 

rH 


CO 

in 

o 


H 

rH 

io 


o 

oo 

CM 

CM 

CO 


CO 

rH 

o 

CM 

05 

e- 

rH 

in 

CO 

05 

rH 

05 

00 

in 

rH 

CD 



O 

05 

CO 

CO 

CM 

rH 

o 

05 

00 

e- 

co 

CD 


r-H 

rH 

rH 


1-1 

rH 

rH 

rH 

rH 








05 

i—l 

05 

CO 

E— 

05 

00 

E- 

CO 

1 — 1 

CO 


CM 


05 

i—l 

CM 

00 

© 

o 


CO 

05 

rji 


05 


i-l 

in 

in 

CM 

CM 

CO 


CM 

CM 

CM 

CM 

CO 



m 

I> 

e~ 

00 




















X> 

CM 


CM 

co 

rH 

CM 

CM 

CM 

in 


oo 

▼H 

05 

oo 

rH 

& 

o 

n 


co 

oo 

in 

rH 


05 

05 

"Cf 

05 

CO 

E- 

in 

CM 

CM 

CO 

co 

05 

03 

eo 

CM 

rH 

o 

o 

00 

oo 

E- 

co 

CD 

in 

CD 

• 


• 

* 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

5 

rH 

rH 

rH 



rH 

rH 

rH 

rH 

rH 








05 

in 

00 

00 

in 


oo 

00 

00 

tH 

E- 

05 

CM 

rH 

E- 

oo 

co 

CO 

CM 


05 

E- 

rH 

00 

o 

03 

oo 

rH 

03 

CO 

E- 

E- 

O 

o 

co 

to 


CM 

rH 

05 

oo 

E- 

CO 

in 


rH 

rH 

rH 

rH 

rH 

rH 

l-H 

H 

rH 








c- 

E- 

E— 

E- 

E- 

00 

E- 

-tf 1 

05 

rH 

CO 

E- 

oo 

o 

tH 

iH 

l—l 

i— 1 

i-H 

00 

o 

O 

03 

O 

O 

rH 

O 

i—l 

CM 

CM 

00 

00 

OO 

oo 

E- 

00 

oo 

E- 

00 

00 

00 

oo 

oo 

oo 

00 


rH 

o 

E- 

CO 

CO 

CO 

CM 


i—l 

in 

o 

CD 

E- 

03 

o 

03 

03 

00 

oo 

05 

03 

00 

05 

oo 

05 

00 

co 

OO 

oo 

E- 

E- 

E- 

03 

03 

05 

03 

03 

05 

03 

05 

03 

03 


988 






TABLE III (Continued) 
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Figure 3. - Critical flow pressure ratio of oxygen, nitrogen, 
argon, and methane computed by isentropic equilibrium 
expansion using corresponding states variables. 
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Figure 4. - Comparison of nitrogen and methane critical 
flow rate data on corresponding states basis. 
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Figure 5. - Comparison of nitrogen and methane critical 
flow pressure ratio data on corresponding states basis. 



